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Abstract. We have conducted deep photometric searches for substellar members of the Praesepe (0.5-1 Gyr) and 
a Orionis (3 Myr) star clusters using the Sloan (' and z' broad-band filters, with the 3.5-m and the 5-m Hale telescopes on 
the Calar Alto and Palomar Observatories, respectively. The total area surveyed was 1177 arcmin 2 and 1122 arcmin 2 towards 
the central regions of Praesepe and <x Orionis, respectively. The 5-cr detection limit of our survey is measured at (' = 24.5 
and z' = 24 mag, which according to state-of-the-art evolutionary models corresponds to masses of 50-55 Mj up (Praesepe) and 
6 M Jup (cr Orionis), i.e., well within the substellar regime. Besides recovering previously known cluster members reported in 
the literature, we have identified new photometric candidates in both clusters whose masses expand the full range covered by 
our study. In cr Orionis, follow-up near-infrared photometry has allowed us to confirm the likely cluster membership of three 
newly discovered planetary-mass objects. The substellar mass function of cr Orionis, which is complete from the star-brown 
\^ | dwarf borderline down to 7 M Jup , keeps rising smoothly with a slope of ar = 0.6*Qj {AN /AM ~ M~ a ). Very interestingly, one 

of the faintest Praesepe candidates for which we have also obtained follow-up near-infrared JH K s photometry nicely fits the 
expected optical and infrared photometric sequence of the cluster. From its colors, we have estimated its spectral type to be 
between L4 and L6. If confirmed as a true Praesepe member, it would become the first L-type brown dwarf (50-60 M Jup ) 
identified in an intermediate-age star cluster. Our derivation of the Praesepe mass function, which is based on state-of-the-art 
evolutionary models, depends strongly on the cluster age. For the youngest possible ages (500-700 Myr), our results suggest 
that there is a deficit of Praesepe brown dwarfs in the central regions of the cluster, while the similarity between the Praesepe 
and cr Orionis mass functions increases qualitatively for models older than 800 Myr. 
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1. Introduction jects to understand the very low-mass stellar and substellar 

populations in the solar neighborhood. 

The knowledge of the substellar mass function of various star- 
forming regions and very young open clusters (< 10 Myr) In this paper, we present the results of deep photometric 
can provide key constrains on our understanding of the for- surveys in two open clusters of different ages: cr Orionis and 
mation mechanisms giving rise to very low-mass stars and Praesepe. Our main objectives are: the identification, for the 
brown dwarfs (masses below 75 Mj up , or 0.072 M Q , Chabrier first time, of the Praesepe unambiguous substellar population, 
& Baraffe 2000J. The study of intermediate-age star clus- and the extension of our previous deep explorations carried out 
ters (100-1000 Myr) is relevant to track down the dynami- in cr Orionis to search for planetary-mass candidates bel ow the 
cal evolution of brown dwarfs and its impact on the shape deuterium burning-mass limit (Zapatero Osorio et al. I2000i. 
of the mass function. It is also important to detail the brown Both clusters have ro ughly solar metallicity (Caballero [2006J 
dwarfs luminosity and effective temperature changes with time. Hambly et al. I1995al > and ar e characterized by v ery litt le in- 
Furthermore, ultracool L- and T-type brown dwarf members of ternal extinction (Lee 119681 Crawford & Barnes 119691 The 
intermediate-age star clusters would become benchmark ob- Praesepe cluster has an age in the interval 0.5- 1 Gyr, and is 

located at a distance of 180 pc (Robichon et al. [T999 ), while 



the youngest cluster, cr Orionis (1-8 Myr, Oliveira et al. 2002; 

Send offprint requests to: Gonzalez-Garcfa, Zapatero Osorio, Zapatero Osorio et al. 120021 with a likely mean age around 
bmgg@laeff.inta.es,mosorio@iac.es 3 Myr) lies at 352 pc (Perryman et al. U997t . 
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Fig. 1. The Palomar (rectangles) and Calar (squares) photometric surveys carried out in the Praesepe (left panel) and 
cr Orionis (right panel) star open clusters. The location of our candidates is indicated with symbols: solid squares correspond to 
candidates with high photometric cluster membership probability, and open circles stand for low probability candidates. Bottom 
images are taken from the Digital Sky Survey (DSS1, red data). [Provided in jpeg format in astro-ph.] 



Previous large scale proper motion studies of Praesepe 
(Hambly et al. U995bl Adams et al. l2002> covered most of the 
lower main sequence of the cluster down to 0.1 M . Deeper 
surveys by Pinfield et al. i 19971 120001 . Magazzu et al. dl998> . 
and Chappelle et al. ( 2005 1 have extended the Praesepe popu- 
lation down to spectral types M9-L0, which corresponds to a 
mass of 0.06-0.08 M Q , i.e., rather close to the substellar border- 
line. Chappelle et al. (2005 1 concluded that Praesepe substellar 
members are likely depleted from the cluster due to dynami- 
cal evolution. Our survey does not cover a significant area of 
Praesepe (thought to be of several tens square degrees), but to 
the best of our knowledge, is the deepest photometric explo- 
ration carried out in the cluster central regions so far. 

We present our optical surveys and cluster member can- 
didate selection criteria in Section 2. Follow-up near-infrared 
photometry (also presented in Section 2) allows us to refine 
photometric membership status of the candidates. In Section 3 
we discuss on the mass of the selected candidates, derive the 
Praesepe and cr Orionis mass functions, and compare them to 
various mass functions available in the literature. We find new 
planetary-mass candidates (6-13 Mj up ) in cr Orionis, and one 
likely Praesepe brown dwarf (50-60 Mj up ) with colors typical 
of L4-L6 spectral types. 

2. Observations 

2.1. Optical photometry 

We have obtained deep CCD images of the young star clus- 
ters Praesepe and cr Orionis using the Sloan i' and z' fil- 
ters and the 5-m Hale Telescope on the Palomar Observatory 
(2002 February 07) and the 3. 5-m telescope on the Calar Alto 
Observatory (2003 December 23, 24). Both telescopes were 
equipped with large format multi-CCD cameras (the Large 
Field Camera, LFC, at Palomar and the Large Area Imager, 
LAICA, at Calar Alto), which cover a substantial area on the 
sky in one single shot. Instrumental details like pixel size, num- 
ber of CCDs, total area covered by each instrument, filter cen- 
tral wavelengths and passbands are summarized in Tabled The 
two filter sets show similar properties. All observations were 
carried out under clear transparency conditions and with see- 
ing values of 1'.'6-T.'7 (Palomar) and 1'.'7-1'.'8 (Calar). 

The log of the optical observations is shown in Table |3 
where we provide the total integration time per filter (i.e., num- 
ber of frames multiplied by individual exposure time) and the 
airmass intervals in which data were collected. Data were ac- 
quired by applying small telescope offsets (dithers of 9" to 
12") from one exposure to the next. This allowed us to cre- 
ate a nightly super-flat-field image as the result of the me- 
dian of all the science frames taken in the same filter. Raw 
data were bias substracted and divided by the corresponding 
normalized superflats. Strong fringing patterns were apparent 



in the z'-band filter; they were corrected down to a level of 
0.5% with the superflats. We aligned all individual flat-fielded 
frames and stacked them altogether to produce the final deep 
i'- and z'-band images free of cosmic rays. Point-like sources 
were identified using the Sextractor software (clearly resolved 
objects were mostly avoided). Some faint sources close to the 
detection limit that were missed by Sextractor were added man- 
ually after visual inspection of the images. 

We depict in Fig.^the regions explored in the Praesepe and 
cr Orionis clusters. There is an overlap of 70 arcmin 2 (Praesepe) 
and 125 arcmin 2 (cr Orionis) between the Palomar and Calar 
surveys. By taking into account these overlapping regions and 
the dither pattern of the observations, the final area covered by 
our optical survey down to the faintest detection limit is 1 177 
arcmin 2 in the Praesepe cluster and 1 122 arcmin 2 in cr Orionis, 
which represents about 2% and 30% of the total cluster areas, 
respectively. The Calar survey imaged regions around the ac- 
cepted centers of mass of the two clusters, while the Praesepe 
Palomar survey was centered close to the location of the pre- 
viously known member Roque Praesepe 1 (RPr 1, Magazzu et 
al. [T998 1, and the Palomar search mapped the Southeast area 
of cr Orionis, in good agreement with the survey of Caballero 
d2006l and Caballero et al. ( 120061 . 

Recently, Holland et al. (2000 1 have claimed that the 
present Praesepe cluster is the result of a recent merging of 
clusters with different ages. The "subcluster" mentioned by 
these authors lies about 3 pc (~1 deg) Southwest from the clus- 
ter centre. We note that our Praesepe survey does not overlap 
with it. Regarding cr Orionis, Jeffries et al. ( 2006 1 have found 
that there are two distinct kinematic groups between the stars 
cr and ( of Orion. Both groups are very young (<10Myr) and 
located at a distance of 350-450 pc; the most numerous group 
is concentrated around cr Orionis, while the population of the 
second component increases towards the North. According to 
Fig. 2 of Jeffries et al. ( I2006> . only the Northwest CCD of the 
Calar survey is located in a region where the number of young 
objects likely belonging to the second group is comparable to 
that of cr Orionis true members. Thus, we do not expect a sig- 
nificant contamination (less than 10%) in our cr Orionis survey. 

The photometric analysis of both the Palomar and Calar 
surveys was performed using routines within the DAOPHOT 
package of IRAF 1 . We have obtained aperture and PSF fitting 
photometry for all sources using circular apertures at least four 
times larger than the average seeing of the frames. Instrumental 
magnitudes have to be converted into observed magnitudes. 
Unfortunately, we did not observe any flux standard star in 
the i' and z' filters in any of the observing campaings. As 



1 IRAF is distributed by National Optical Astronomy 
Observatories, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under contract to the National 
Science Foundation, USA 
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Fig. 2. Optical color-magnitude diagramas of our Palomar (left) and Calar (right) surveys in the Praesepe star cluster. Selected 
Praesepe candidates are plotted as filled circles. Previously known members, which are labeled, define the cluster photometric 
sequence depicted with a dashed line. The solid lines indicate the limiting magnitudes of our surveys; completeness is about 1 
mag brighter. Masses for an age of 500 Myr are indicated to the right side of the diagrams. Note that only PSF photometric errors 
are plotted. 
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Fig. 3. Optical color-magnitude diagramas of our Palomar (left) and Calar (right) surveys in the <x Orionis star cluster. Selected 
<x Orionis candidates are plotted as filled circles. Previously known members, which are labeled, define the cluster photometric 
sequence depicted with a dashed line. The solid lines indicate the limiting magnitudes of our surveys; completeness is about 1 
mag brighter. Masses for an age of 3 Myr are indicated to the right side of the diagrams. Note that only PSF photometric errors 
are plotted. 
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Table 1. Optical instruments employed in our survey. 















Sloan 


(' filter 


Sloan z' filter 


Telescope 


Observatory 


Camera 


Detectors 


Pixel size 


Area 


A c 


AA 


A c AA 










C'/pix) 


(arcmin 2 ) 


(nm) 


(nm) 


(nm) (nm) 


5-m Hale 


Palomar 


LFC a 


6x(1024x2048) 


0.360 


450 


768 


154 


900 180 


3.5-m CAHA 


Calar Alto 


LAICA 


4x (4096x4096) 


0.225 


940 


770 


150 


920 160 



a The original size of each CCD detector is 2048x4096, and the pixel size is 07 18. We applied a 2x2 binning to all LFC CCDs. 
Table 2. Log of optical observations. 



Cluster 


Survey 


Obs. date 


Filter 


Exp. 




Airmass 






(UT) 




(s) 


(mag) 




cr Orionis 


Palomar 


2002 Feb 07 


i' 


46x200 


23.2, 24.5 


1.38-1.39 




Palomar 


2002 Feb 07 


z' 


21x200 


22.5, 24.0 


1.40-1.93 




Calar 


2003 Dec 23 


i' 


26x500 


23.8, 24.5 


1.30-1.93 




Calar 


2003 Dec 24 


z' 


26x600 


22.8, 24.0 


1.30-1.95 


Praesepe 


Palomar 


2002 Feb 07 


V 


28x200 


23.8, 24.8 


1.03-1.23 




Palomar 


2002 Feb 07 


z' 


21x200 


23.3, 24.3 


1.25-1.80 




Calar 


2003 Dec 23 


i' 


9x500 


23.3, 24.5 


1.05-1.13 




Calar 


2003 Dec 23-24 


z' 


18x500 


22.8, 24.0 


1.05-1.44 



Table 3. Previously known Praesepe candidates in our survey. 



Object 


z' 


(' - z' 


r 


Reference 




WFC210 


17.58±0.20 


1.13+0.28 


18.02 


Chappelle et al 


. (20051 


PHJ20 


17.65±0.20 


1.21+0.28 


18.48 


Hodgkin et al. I 


T3391 


WFC193 


17.81±0.20 


1.13+0.28 


18.21 


Chappelle et al 


. <2005l 


PHJ18 


18.43+0.20 


1.64+0.28 


19.63 


Hodgkin et al. 1 1999 ( 


WFC88 


19.06±0.20 


1.51+0.29 


19.83 


Chappelle et al 


. <2005l 


RPrl 


19.71±0.20 


1.53+0.29 


21.01 


Magazzii et al . 


(T9981 


WFC60 


19.86±0.20 


1.66+0.28 


20.79 


Chappelle et al 


. <2005l 


WFC24 


20.32±0.20 


1.66+0.29 


21.19 


Chappelle et al 


. <2005l 



a /-band magnitude from the literature. 



on January 2006, no public data overlapping the regions of 
our survey were available to us from the Sloan Digital Sky 
Survey. Hence, we calibrated instrumental photometric data us- 
ing previously identified Praesepe and cr Orionis members that 
also appear in the Palomar and Calar surveys. For these ob- 
jects, Magazzu et al. ( 119981 1. Hodgkin et al. ( 1999), Zapatero 
Osorio et al. ( 2000} , Bejar et al. J2DM . Chappelle et al. ( I20031 . 
Kenyon et al. J2005> . and Caballero et al. i2004l 120061 pro- 
vide /-band photometry, and for many the (R - /) color is also 
available. Tables|3]and0]show the complete list of known clus- 
ter members in our survey along with their /-band magnitudes 
taken from the literature and from the very recent calibration 
of Bejar et al. (2006). We note that Bejar et al. (2006) provide 
a consistent and homogeneous optical photometric calibration 
for a large area of the cr Orionis cluster. 

The known objects of Tables[3]and|4]are all red with likely 
spectral types cooler than mid-M. Their /-band magnitudes 
were transformed into the Sloan i' and z' data by applying the 
following equations: 

i' = / + 0.398 + 0.143(7?-/) (1) 



j"-z' = 1.22(i?-7)-1.47 (2) 

which are valid for red objects with (R- 1) > 0.75 mag (Eq.[Q 
and (R - I) > 1.65 mag (Eq.0. Smith et al. d2"0"02l provide 
the Sloan photometry of a large sample of equatorial standard 
stars drawn from the Landolt ( 1992 1 catalog. We have related 
their i' data to the Landolt R and / magnitudes to obtain eq. ^ 
Equation|5]was derived from the relations of Smith et al. (2002 1 
and the theoretical calculations by Fukugita et al. {I996l. The 
i' and z' values computed for the known cluster members were 
then compared to the objects' instrumental magnitudes from 
our search, and used to define the zero point calibration of the 
data. We note that because of this procedure, which strongly 
relies on the colors of the calibrating sources, the i' and z' mag- 
nitudes of objects with blue colors are not reliable at all in our 
work. However, our main interest is the study of the very red 
population of the Praesepe and <x Orionis clusters, for which 
the calibration should be reasonable. The ex Orionis Palomar 
survey was calibrated using objects from the overlapping re- 
gion with the Calar survey. The i'- and z'-band dispersion of 
the adopted photometric calibration for the two clusters turns 
out to be +0.2 mag, and the error associated to the zero points 



B. M. Gonzalez-Garcfa et al.: A search for substellar members in the Praesepe and a Orionis clusters 
Table 4. Previously known cr Orionis candidates in our survey. 
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Object 




z' 


i' — z' 


P 


Reference 


SOri 13 




15.66±0.21 


1.05+0.30 


16.05 


Bejaret al. 1 2006 1 


S Ori 14 




15.92±0.21 


1.26±0.29 


16.35 


Bejar et al. 119991 


S Ori 19 




16.11+0.21 


1.08+0.29 


16.62 


Bejar et al. <2006l 


S Ori 23 




16.43±0.21 


1.14+0.29 


16.88 


Bejar et al. <2006> 


SOri 24 




16.48±0.20 


1.10+0.29 


16.90 


Bejaret al. 1 2006 1 


S Ori 28 




16.74±0.21 


1.24+0.29 


17.23 


Bejaret al. 1 2006 1 


SOriJ053922.2- 


-024552 


16.78±0.21 


1.08+0.29 


17.05 


Caballero et al. <2004t 


S Ori 32 




16.89±0.21 


1.26+0.29 


17.40 


Bejar et al. <2006> 


KNJ2005 72 




17.01+0.21 


1.31+0.29 


17.82 


Kenyon et al . 120051 


KNJ2005 70 




17.13±0.21 


1.29+0.29 


17.81 


Kenyon et al . |2005 I 


KNJ2005 71 




17.26±0.21 


1.12+0.29 


17.82 


Kenyon et al . (2005 l 


SOri 38 




17.34±0.24 


1.32+0.35 


17.70 


Bejar et al. <2006> 


SOriJ053926.8- 


-022614 


17.69±0.21 


1.15+0.30 


18.25 


Bejar et al. <2006> 


S Ori J0539 18.1 ■ 


-025257 b 


17.99±0.21 


1.66+0.29 


18.64 


Caballero et al. <2004t 


SOriJ053948.1- 


-022914 


18.06±0.21 


1.20+0.30 


18.52 


Bejar et al. <2006> 


S Ori 42 




18.33+0.21 


1.70+0.30 


19.04 


Caballero et al. <2004t 


SOri 44 




18.94±0.21 


1.29+0.29 


19.39 


Bejar et al . 119991 


SOri 45 




18.99±0.21 


1.58+0.30 


19.49 


Caballero et al. <2004t 


SOriJ053929.4- 


-024636 


19.12±0.21 


1.68+0.29 


19.73 


Caballero et al. i2004l 


SOriJ053946.5- 


-022423 


19.09±0.21 


1.59+0.29 


20.14 


Bejar et al. <2"0oTl 


S Ori 52 




19.76±0.21 


1.78+0.29 


20.55 


Bejar et al. <2006l 


SOri 53 




20.42+0.22 


1.85+0.31 


20.96 


Caballero et al. <2004l 


SOriJ053956.8- 


-025315 


20.49±0.21 


1.86+0.30 


21.25 


Caballero et al. <2Ti04l 


SOriJ053858.6- 


-025228 


21.28±0.24 


1.98+0.35 


22.17 


Caballero et al. iSoTBl 


S Ori 67 




22.05±0.22 


1.78+0.35 


23.04 


Bejaret al. 1 2006 1 



a /-band magnitude from the literature. 

b It is detected at both Palomar and Calar surveys. We provide the mean i' and z' photometry. The i' magnitudes differ by 0.3 mag suggesting 
photometric variability and supporting the conclusions on the source variable nature reported by Caballero et al. 1 2004 1. 



Table 5. New Praesepe candidates from our survey 



IAU name 


RA (J2000) 

^li m s j 


DEC (J2000) 

C ' ") 


Palomar 

z' 


survey 
(' - z' 


Calar survey 
Zf i' - 


-z 1 


P 


Phot. Prob. 


PraeJ083850.6+192317 


08 38 50.55 


+19 23 16.5 






17.57+0.20 


1.08d 


=0.28 


18.20 


high 


PraeJ084051.2+191010 


08 40 51.24 


+ 19 10 09.5 


17.75+0.20 


0.93+0.28 








18.41 


low 


PraeJ084108.0+191401 


08 41 08.00 


+ 19 14 00.6 


18.42+0.20 


1.04+0.28 








19.20 


low 


PraeJ083905. 4+ 192032 


08 39 05.38 


+ 19 20 31.5 






19.16+0.20 


1.26=: 


=0.28 


20.06 


low 


PraeJ084 130.4+ 190449 


08 41 30.43 


+ 19 04 48.5 


19.21+0.20 


1.16+0.28 








20.12 


low 


PraeJ0841 15.8+191743 


08 41 15.79 


+ 19 17 43.0 






19.39+0.20 


1.19= 


=0.28 


20.33 


low 


PraeJ084114.5+191110 


08 41 14.46 


+ 19 11 10.1 


19.42+0.20 


1.14+0.28 








20.36 


low 


PraeJ084126.5+195200 


08 41 26.54 


+ 19 51 59.7 






19.78+0.20 


1.38= 


=0.28 


20.78 


low 


PraeJ084108.1 + 191111 


08 41 08.13 


+ 19 11 10.8 


20.07+0.20 


1.42+0.29 








21.12 


low 


PraeJ084039.3+ 192840 


08 40 39.28 


+ 19 28 39.6 






20.26+0.20 


1.78= 


=0.28 


21.34 


high 


PraeJ084052.5+ 190504 


08 40 52.52 


+ 19 05 03.5 


20.73+0.20 


1.42+0.29 








21.89 


low 


PraeJ084122.6+192419 


08 41 22.58 


+ 19 24 19.0 


21.33+0.25 


1.55+0.35 


21.21+0.20 


1.61d 


=0.28 


22.52 


low 


PraeJ084049.1 + 195913 


08 40 49.14 


+ 19 59 12.6 






21.30+0.20 


1.71d 


=0.29 


22.55 


low 


PraeJ084044.9+200025 


08 40 44.90 


+20 00 25.3 






21.58+0.20 


1.64= 


=0.29 


22.88 


low 


PraeJ084045.2+ 192532 


08 40 45.24 


+ 19 25 32.3 






21.72+0.20 


1.68= 


=0.29 


23.05 


low 


PraeJ083910.2+195610 


08 39 10.20 


+ 19 56 10.1 






21.86+0.20 


1.65= 


=0.29 


23.21 


low 


PraeJ084111.9+192043 


08 41 11.90 


+ 19 20 43.3 


22.33+0.22 


1.79+0.32 


22.33+0.20 


2.27= 


=0.43 


23.76 


high 


PraeJ083904.7+ 194959 


08 39 04.72 


+ 19 49 58.6 






22.36+0.20 


1.87= 


=0.31 


23.80 


low 


PraeJ084057.3+ 190258 


08 40 57.37 


+ 19 02 57.7 


22.49+0.27 


1.84+0.35 








23.95 


low 


PraeJ0841 11.2+191349 


0841 11.15 


+ 19 13 49.1 


22.75+0.25 


1.82+0.36 








24.25 


low 



a Estimated from the z' photometry. 
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Table 6. New cr Orionis candidates. 



IAU name 




RA (J2000) 


DEC (J2000) 

C ' ") 


Palomar 

z ' 


survey 
i' - z' 


Calar survey 
z' i' - z' 


r 


Phot. Prob. 


SOriJ053724.5- 


-021856 b 


05 37 24.46 


-02 18 55.7 






16 58+0 21 


1.12±0.29 


17.07 


high 


SOriJ053855.0- 


-024034 


05 38 54.95 


-02 40 33.6 


17.69+0.21 


1.64±0.29 






18.24 


high 


SOriJ053734.5- 


-025559 


05 37 34.51 


-02 55 58.5 






17.79±0.21 


1.35±0.30 


18.35 


high 


SOriJ053908.4- 


-023902 


05 39 08.38 


-02 39 01.9 


18.03±0.22 


1.22+0.31 






18.60 


low 


SOriJ053912.6- 


-025724 


05 39 12.60 


-02 57 23.9 


18.69+0.21 


1.30±0.29 






19.29 


low 


SOriJ053731.9- 


-021634" 


05 37 31.89 


-02 16 33.8 






19.03±0.21 


1.49±0.29 


19.65 


high 


SOriJ053946.3- 


-022631 


05 39 46.33 


-02 26 31.3 






19.70±0.29 


1.5 1+0.39 


20.36 


low 


SOriJ053913.9- 


-021621 b 


05 39 13.94 


-02 16 21.2 






19.80±0.21 


1.73±0.29 


20.47 


high 


S Ori J0538 13.5- 


-021934 b 


05 38 13.54 


-02 19 34.2 






20.13±0.21 


1.48±0.29 


20.83 


low 


SOriJ053758.0- 


-025146 


05 37 57.98 


-02 51 46.1 






20.63±0.26 


1.59+0.38 


21.38 


low 


SOriJ053957.4- 


-025006 


05 39 57.40 


-02 50 06.2 






21.26±0.22 


1.95+0.32 


22.10 


high 


SOriJ053919.2- 


-025910 


05 39 19.23 


-02 59 10.4 


21.60+0.28 


1.84±0.38 






22.50 


low 


SOriJ054007.8- 


-022234 


05 40 07.75 


-02 22 34.3 






21.65±0.22 


1.79±0.33 


22.56 


low 


SOriJ053844.5- 


-0255 14 c 


05 38 44.54 


-02 55 14.3 


21.66±0.27 


2.02±0.38 






22.57 


high 


SOriJ054008.5- 


-02455 l c 


05 40 08.51 


-02 45 50.6 






21.85±0.22 


1.91 ±0.44 


22.80 


low 


SOriJ053932.4- 


-025220 c 


05 39 32.43 


-02 52 20.2 


21.99±0.22 


1.85±0.40 


2 1.91 ±0.22 


2.13±0.40 


22.87 


high 


SOriJ053732.3- 


-022842 


05 37 32.32 


-02 28 42.1 






22.09±0.23 


1.98±0.37 


23.09 


low 


S Ori J0538 16.3- 


-022727 


05 38 16.38 


-02 27 27.1 






22.18±0.24 


1.85±0.37 


23.20 


low 



a Estimated from the z' photometry. 
b Also in Bejar et al. <2006> . 
c Also in Caballero et al. (2006 1. 



is 0.06 mag. The consistency of the calibration between the 
Palomar and Calar surveys in Praesepe was checked by exam- 
ining the sources in the overlapping area: for any of the filters, 
we found no obvious offset larger than the uncertainty of the 
calibration. The relatively large dispersion in the obtention of 
the zero points is probably due to the use of various photomet- 
ric sources for the /-band magnitudes, to the dispersion intro- 
duced by the transformation eqs. [2 and [2 and to the possible 
intrinsic variability of some objects. 

We provide in Tables [3] (Praesepe) and 0] (cr Orionis) the 
final i' and z' magnitudes of the cluster members in common 
with the literature as obtained from our survey. Photometric 
error bars account for both the PSF error provided by the 
IPsAF routines and the uncertainty in the photometric calibra- 
tion. Limiting (4-5-cr) and completeness (10-cr) magnitudes 
for each filter and survey are provided in Tabled In short, our 
search may be sensitive to objects down to i' = 24.5 and z' = 24 
mag, while the survey completeness is determined at i' ~23.5 
and z' ~ 23 mag. However, because we impose that the cluster 
candidates are detected in both filters and the faintest cluster 
members should have very red colors (/' - z' ~ 2), our survey 
is actually limited by the blue /'-band. 

2.2. Identification of optical candidates 

The color-magnitude diagrams of our optical survey are illus- 
trated in Figs. |2 (Praesepe) and[3](cr Orionis). We identify the 
previously known cluster members with labels in the Figures. 
These objects are clearly separated from the field sources in 
the diagrams, and define a photometric sequence character- 
istic of each cluster. The extrapolation of the z' vs. (/' - z') 



cluster sequence down to the limit of our survey is plotted as 
dashed lines. As can be seen from the Figures, the slope of the 
Praesepe and cr Orionis dashed lines is quite similar provided 
the fact that true cluster members become increasingly redder 
for fainter magnitudes. The limiting magnitudes of both filters 
are also indicated by solid lines in the Figures. 

To select cluster member candidates we have plotted an 
imaginery straight line that is parallel to the cluster photometric 
sequence and shifted 0.4 mag towards the blue. For bluer col- 
ors the number of detections increases very drastically, proba- 
bly indicating that these are field objects. Additionally, this line 
coincides with the apparent lower envelope delineated by the 
previously known candidates of the cr Orionis cluster (e.g., see 
the location of S Ori 44 and S Ori 67 in Fig. The photomet- 
ric dispersion observed in other optical colors among low-mass 
proper motion members of the Pleiades and Praesepe clusters 
is about 0.6-0.7 mag. 

A total of 28 candidates have been extracted from our sur- 
vey of the Praesepe cluster, 20 of which are new. Regarding 
cr Orionis, 43 candidates have been found to match our selec- 
tion criteria, 18 of which have no previous mention in the liter- 
ature. Tables[5]and[6]provide IAU names, coordinates and the 
optical photometry of all the new candidates. The astrometric 
plate solution has been obtained for all the CCDs by compar- 
ison of the observations with the USNO B-l catalog (Monet 
et al. 120031 1. Possible spatial distortions have been modelled 
in a standard astrometric reduction procedure using the task 
PLTSOL of IRAF. The quadratic mean (rms) of the deviations 
between the catalog and the computed coordinates for the trans- 
formation stars is +0.6" (right ascension) and +0.5" (declina- 
tion) for both the Palomar and Calar surveys. 
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Table 7. Near-infrared photometry of a few new Praesepe candidates. 



Object 


Teles. 


J 




H 


K, 


z'-J 


J - K s 


Phot. prob. a 


PraeJ083850.6+192317 b 


2MASS 


16.374 


:0.11 


16.01+0.18 


15.40+0.19 


1.20+0.23 


0.97+0.12 


high / high 


PraeJ084051.2+191010 c 


2MASS 


16.75+ 


:0.15 


15.76+0.15 


15.59+0.23 


1.00+0.25 


1.16+0.27 


low / high 


WFC24 d 


TCS 


17.71i 


:0.33 






2.6+0.4 




high / high 


PraeJ0841 11.9+192043 


WHT 


19.85+ 


:0.07 


18.95+0.11 


18.28+0.10 


2.48+0.21 


1.57+0.12 


high / high 



a Optical / near-infrared photometric cluster membership probability. 

b 2MASS J08385059+1923164. 

c 2MASS J08405122+1910106. 

d J= 17.99±0.03, H= 17.30+0.02, K = 16.78+0.02 (Chappelle et al. 2005 l. 



Table 8. Near-infrared photometry of new <x Orionis candidates. 



Object 




J 


z'-J 


Phot. prob. a 


Ref. 


SOriJ053724.5- 


-021856 


15.36+0.09 


1.22+0.23 


high / high 


Bejar et al. ( 2006 1 


SOriJ053734.5- 


-025559 b 


15.71+0.06 


2.08+0.22 


high / high 


2MASS 


SOriJ053855.0- 


-024034° 


15.92+0.07 


1.77+0.22 


high / high 


2MASS 


SOriJ053908.4- 


-023902 d 


16.46+0.10 


1.57+0.23 


low / high 


2MASS 


SOriJ053731.9- 


-021634 


17.46+0.08 


1.57+0.23 


high / low 


Bejar et al. i2006l 


SOriJ053913.9- 


-021621 


17.57+0.10 


2.23+0.24 


high / high 


Bejar et al. t2006> 


SOriJ053813.5- 


-021934 


18.56+0.20 


1.57+0.28 


low / low 


Bejar et al. i2006l 


SOriJ053844.5- 


-025514 


19.31+0.08 


2.35+0.28 


high / high 


Caballero <2006l 


SOriJ054008.5- 


-024551 


19.43+0.16 


2.43+0.27 


low / high 


Caballero et al. <2006l 


SOriJ053932.4- 


-025220 


19.54+0.08 


2.37+0.24 


high / high 


Caballero <2006l 


SOriJ053919.2- 


-025910 


19.83+0.09 


1.77+0.29 


low / low 


Caballero et al. 1 2006 1 



a Optical / near-infrared photometric cluster membership probability. 

b 2MASS J05373457-0255589, H = 14.86+0.06, K s = 14.39+0.08 mag. 

c 2MASS J05385492-0240338, H = 15.17+0.06, K s = 14.71+0.11 mag. 

d 2MASS J05390846-0239013, H = 15.72+0.10, K s = 15.19+0.13 mag. 



Additionally, because our selection criterion may be quite 
generous, for each new candidate we have asigned a flag on 
the cluster membership probability column of Tables |5] and |6] 
based solely on the object's proximity to the expected optical 
photometric sequence. If the (/' - z') color of the candidate lies 
within 0.2 mag from the expected index given by the dashed 
lines of the color-magnitude planes of Figs. |2 and then the 
flag turns to "high" (indicating high photometric cluster mem- 
bership probability), otherwise it is "low". In any case, further 
follow-up data are required to confirm the true membership of 
all the candidates in their respective clusters. The eighth col- 
umn of Tables|5]and[6]provide the broad-band / magnitudes of 
the new objects, which have been estimated from their Sloan z' 
photometry. Finder charts using the z' frames for both "high" 
and "low" probable candidates are provided in Figs. l9ltol 1 21 

A few new Praesepe and <x Orionis candidates have been 
found in the overlapping regions of the surveys. In gen- 
eral, the Calar and Palomar photometry do agree well, ex- 
cept for the faintest candidates (S Ori J053932.4-025220 and 
PraeJ084111. 9+192043), where large discrepancies (up to 
0.45 mag) can be observed in the /'-band data; on the contrary, 
the z' magnitudes are quite alike within less than 1-cr the er- 
rors. This is explained by the fact that these objects appear to 
be notably red and have /'-band magnitudes indeed very close 



to the detection limit of our survey, where we expect large pho- 
tometric uncertainties. 

2.3. Near-infrared photometry 

A few of the brightest new Praesepe and cr Orionis candi- 
dates of Tables[5]and[6]were detected by the 2 Micron All-Sky 
Survey (2MASS, Skrutskie et al. 12006V Their JHK S magni- 
tudes and associated error bars are given in Tables and [8] 
The 2MASS photometry shows rather large uncertainty be- 
cause these objects are close to the detection limit of the all- 
sky survey. Nevertheless, both optical and near-infrared data 
support them as red sources. 

Additionally, we carried out point observations of two 
of our optical Praesepe candidates using the J, H, and 
K s filters and the 1.5-m Carlos Sanchez Telescope (TCS, 
Teide Observatory) and the 4.2-m William Herschel Telescope 
(WHT, Roque de los Muchachos Observatory). The date of the 
observations were as follows: 2005 March 07 at the TCS, and 
2006 March 20 and 22 at the WHT. We used the HgCdTe, 
l"-pixel detector at the Cassegrain focus of the TCS, and 
the LIRIS camera (equipped with a 0.25"-pixel, 1024x1024 
HAWAII detector) at the WHT. Observations were performed 
in the broad-band J, H, and K s filters. The total integration time 
was 2400 s for WFC24 (/-band), 640 s (J) and 800 s (H, K s ) 
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i'-z' i'-z' 

Fig. 4. Our Praesepe (left panel) and cr Orionis (right panel) candidates are plotted in the optical color-magnitude diagram. Filled 
circles stand for the newly discovered cluster candidates with high membership probability, while previously known cluster 
members are plotted as open circles. Our "less" probable member candidates are depicted with open triangles (they lie in a 
sequence that appears slightly bluer than those of the clusters). Overplotted onto the data are the 500-Myr (Praesepe) and 3-Myr 
(cr Orionis) isochrones (solid lines) from the Lyon group (Chabrier et al. 2000; Baraffe et al. 1998 2003 1. Masses in Jovian units 
are indicated to the right side of the diagrams. The field sequence of late-M and L-type dwarfs at the distance of the Praesepe 
cluster is plotted as a dotted line in the left panel (spectral types are labeled). Note that only PSF photometric errors are shown. 



for PraeJ0841 1 1.9+192043, where each final image actually 
consisted of many (32-240) short individual exposures of 5- 
20 s. Raw near-infrared data were reduced in a standard way 
for these wavelengths, which include sky substraction and fiat- 
fielding. Weather conditions during the campaigns were always 
photometric at the WHT, and some cirrus were affecting the 
TCS observations. The photometric analysis of the two candi- 
dates was carried out using routines within IRAF. Because of 
the faintness of the objects, we obtained PSF fitting photom- 
etry. Instrumental magnitudes were transformed into observed 
magnitudes using 2MASS sources that were present within the 
field of view of the detectors, i.e., the targets and the calibra- 
tors were observed simultaneously. Our final meaurements are 
given in Table along with their associated error bars, which 
take into account the uncertainties of the PSF photometry and 
the photometric calibration. We note that WFC 24 was detected 
quite close to the detection limit of the TCS near-infrared im- 
age. Nevertheless, our /-band measurement is in quite well 
agreement with the data by Chappelle et al. ( 2005 ) within the 
expected errors. 

Our second near-infrared target in the Praesepe cluster is 
Prae J0841 1 1 .9+192043, one of the faintest and most attractive 
photometric Praesepe candidates found in the optical survey. 
The near-infrared data do confirm it as a very red, unresolved 
object (average seeing of the WHT observations is 0785). This 
object will be discussed further in the following Section. 



Near-infrared photometry of eight new cr Orionis candi- 
dates are reported by Caballero ( 2006 ), Caballero et al. ( 2006 1, 
and Bejar et al. (2006 1, where the details of data obtention and 
reduction can be found. We summarize in Table[8]their /-band 
magnitudes and associated uncertainties. As deduced from the 
(z' - /) colors, all objects turn out to be red, but not all of 
them are red enough to fit the cr Orionis near-infrared sequence 
of true cluster members. Similarly to what we did using the 
optical data, we have also assigned a "low" or "high" near- 
infrared photometric cluster membership probability flag for 
each one of the Praesepe and cr Orionis candidates in Tables 
and Cluster membership based on the various optical and 
near-infrared color-magnitude diagrams will be discussed next. 

3. Discussion 

3.1. Cluster membership and mass of the candidates 

Our photometric survey was designed to detect in both i' and 
z' filters Praesepe brown dwarfs down to about 50 times the 
mass of the planet Jupiter (1 M = 1047 M Jup ), or 55 M Jup 
(completeness), i.e., well below the substellar mass limit de- 
fined at 75 Mj up (Chabrier & Baraffe I20001 This prediction 
is based on the state-of-the-art evolutionary models by Baraffe 
et al. (1998 1, an age of 500 Myr, and a distance modulus of 
(m - M) = 6.28 mag (180 pc, Robichon et al. 1999}. For the 
older age of 1 Gyr, the mass limit of our survey is estimated 
at 60 Mj up , still in the substellar regime. In terms of spectral 



B. M. Gonzalez-Garcfa et al.: A search for substellar members in the Praesepe and a Orionis clusters 



9 



CD 



co 



o 

03 



-i 1 1 1 1 1 1 1 1 1 1 1 1 1- 

Praesepe 



o o 
o o 

O o 

o o 



o u o 




M/M 

75 
63 

52 
42 



Jup 



_l I I I I I I I I I I I l_ 



1.5 



J - K s 

Fig. 5. Near-infrared color-magnitude diagram of Praesepe 
members. Our "low" and "high" new optical candidates with 
JK S photometry are plotted as open triangles and filled circles, 
respectively, along with their associated error bars. Previously 
known members and data from the literature are shown with 
open circles. The field sequence at the distance of the cluster is 
plotted as a solid line (spectral types are labeled). Masses ac- 
cording to an age of 500 Myr and the models by Chabrier et al. 



2000 1 are indicated to the right side of the diagram. 



type, these brown dwarfs would show mid- to late-L spectral 
appearance with surface temperatures of about 1500 K (Vrba 
et al. 120041 . Similarly, our <x Orionis survey is sensitive to 
the same spectral interval, which corresponds to masses of 6- 
7 Mj up (limit) and 8 Mj up (completeness) at the age (3 Myr) and 
distance (352 pc) of the cluster. 

There are various new Praesepe and <x Orionis candidates 
discovered in our survey that have been labeled with a "high" 
flag in the optical photometric cluster membership columns 
of Tables [5] and [6] All of them are plotted as filled circles in 
Fig.0] where previously known cluster members from the lit- 
erature are also shown. Overplotted onto the data are the 500- 
Myr (Praesepe) and 3-Myr (cr Orionis) isochrones obtained 
from the NextGen (3000-2000 K, Baraffe et al. [1998), dusty 
(2000-1300 K, Chabrier et al. 2000) and COND (1300-800 K, 
Baraffe et al. 2003 1 evolutionary models. To derive the pre- 
dicted i' and z' magnitudes from the theoretical r e ff's and lu- 
minosities, we have used the relations between the observed 
Sloan (/' — z') colors of field late-M and L dwarfs and their mea- 
sured surface temperatures and /-band bolometric corrections 
compiled from the works by Knapp et al. (2004 1, Vrba et al. 
(20041, and Dahn et al. (2002J. As can be seen from Fig.g] the 
models nicely reproduce the shape and trend of the Praesepe 
and cr Orionis optical photometric sequences. Given the age of 
Praesepe (0.5-1 Gyr), which is close to that of the field, the 
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Fig. 6. Optical and near-infrared color-magnitude diagram of 
cr Orionis members. Our candidates with i'z' high photometric 
membership probability are plotted as filled circles, while open 
triangles stand for low probability ones. Previously identified 
cluster members are shown with open circles. Overplotted onto 
the data is the 3-Myr isochrone from Chabrier et al. (2000) and 
Baraffe et al. (2003). Masses are labeled to the right. 



cluster sequence also matches within the expected photometric 
errors the sequence delineated by field very low-mass stars and 
brown dwarfs as it is shown in the Figure. The theoretical 500- 
Myr isochrone appears slightly overluminous as expected for a 
younger age. The largest differences in luminosity are observed 
at colors (/' — z') > 3 mag. 

Figures[5]and[6]illustrate the near-infrared color-magnitude 
diagram of the Praesepe cluster and the optical-near infrared 
diagram of cr Orionis, respectively. As before, the field se- 
quence moved to the distance of Praesepe is overplotted in 
Fig. |3 It provides a good guideline of the expected pho- 
tometric trend that true Praesepe members should follow at 
J and K s magnitudes. One of the new brightest candidates, 
Prae J083850.6+192317, nicely fits the cluster photometric 
sequence at both optical and near-infrared wavelengths (see 
Figs.HandHl, as also does Prae J08405 1.2+ 19 1010, although 
this candidate was labeled with a "low" flag in Table [5] These 
objects are likely Praesepe stars with masses estimated at 0.10- 
0.15 M e . 

Another probable photometric Praesepe candidate is 
Prae J084039. 3+ 192840, which has similar brightness and col- 
ors than those of the faintest cluster members identified by 
Chappelle et al. (2005). Because of its similarity to RPr 1 



(Magazzu et al. [T998 ) and based on its optical color, we es- 
timate its spectral type to be between M8 and LI, i.e., in 
the M-L transition. According to the evolutionary models, 
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Fig. 7. Images (2' x 2' in extent) of PraeJ0841 1 1.9+192043 taken at four different wavelengths. [Provided in jpeg format in 
astro-ph.] 



Prae J084039. 3+ 192840 sits on top of the expected star-brown 
dwarf frontier of the Praesepe cluster. 

Very interestingly, we have identified an outstanding 
Praesepe candidate member that is signficanfly fainter and red- 
der than any other candidate: Prae J0841 1 1.9+192043, which 
is detected as a very red source in both the Palomar and Calar 
optical surveys. As can be seen from Fig.|4j its optical colors 
are consistent with mid-L spectral type, and models predict a 
mass of 50-60 Mj up for it (0.5-1 Gyr). The near-infrared pho- 
tometry of Table and the compilation of images depicted in 
Fig0confirm its cool nature. Prae J0841 1 1 .9+192043 does fol- 
low the expected Praesepe cluster sequence both at optical and 
near-infrared wavelengths (Figs.|4]and|5|i. We can estimate the 
number of L4-L6 field dwarfs that may appear in our survey 
with similar magnitudes to those of Prae J0841 1 1.9+192043 
by taking into account the local number density of these ob- 
jects in the solar neighbourhood. Cruz et al. (|2003) found an 
L0-L4 dwarf space density of 1.9xl0~ 3 objects pc~ 3 , while for 
T dwarfs (which are cooler than the L-type objects), Burgasser 
(2001 1 derived a bias-corrected value of 21xl0~ 3 objects pc -3 . 
In qualitative terms, there are more T dwarfs than L dwarfs. We 
will adopt a mean density of 10xl0~ 3 dwarfs with mid-L types 
per cubic parsec, which is intermediate between the early-L 
and T classes. For the area explored in the Praesepe survey, 
we would have expected 0.4 such field objects contaminating 
our optical search, i.e., less than 1. There is a ~70% proba- 
bility that Prae J084 11 1.9+192043 is a true Praesepe member. 
We also note that there is a negligible probability for this ob- 
ject being a background very red galaxy. On the one hand, it 
remains unresolved in our good-seeing WHT near-infrared im- 
ages. On the other hand, Fig. 10 of Metcalfe et al. ( 2006 1, which 
displays the optical near-infrared color-magnitude diagram of 
galaxies in the interval H - 16-22.5 mag, shows an empty area 
around the location defined by the colors and magnitudes of 
Prae J0841 1 1.9+192043. Nevertheless, further data are needed 
to unambiguously confirm its cluster membership, for example, 
lithium detection in its optical spectrum and proper motion. 

Regarding the cr Orionis cluster, we note that with a few 
exceptions, there is a one to one relation in the optical and 
near-infrared photometric cluster membership probability flags 
given in Table [8] By considering the optical and near-infrared 
photometric data of both known and new cr Orionis candidates 
(36 objects out of a total of 43 candidates have photometry at 
various wavelengths), we estimate that the success rate on the 
identification of likely cluster members in our survey is about 
93% among the "high" candidates, and about 38% among the 
"low" ones. Hence, to build reliable cluster luminosity and 
mass functions we discard the "low" optical candidates un- 
less they are proved to follow the near-infrared sequence of 
cr Orionis. 

Of special interest is the population number in the 
planetary-mass regime of the cr Orionis cluster. In our sur- 
vey, we have recovered two previously known cluster candi- 
dates whose masses could be below the deuterium burning- 



mass limit (z' > 20.5 mag, J > 18 mag). Additionally, we 
have identified three high-probable and six low-probable new 
planetary-mass candidates. Near-infrared photometry is avail- 
able for four of them (two "high" and two "low"). As can 
be seen from Fig. [6] one low-probable planetary-mass can- 
didate (SOriJ053919.2-025910) does not show an infrared 
color consistent with membership in cr Orionis, and thus, it 
is not considered further in this paper. On the other hand, 
SOriJ053844.5-025514and SOri J053932.4-025220, which 
were given a "high" flag based on their optical data (see the 
right panel of Fig.0}, also fit the near-infrared cluster sequence. 
This suggests that they are likely planetary-mass members of 
the cr Orionis cluster. In addition, S OriJ054008. 5 -024551 
shows infrared colors expected for true cluster members de- 
spite the fact it was originally classified as a low-probable op- 
tical candidate. 

According to modern evolutionary models, 
SOriJ053844.5-025514, S Ori J053932.4-025220, and 
S OriJ054008. 5 -024551 may have masses between 6 and 
10 Mj up if their cluster membership is finally confirmed. In 
total, there are 6 cr Orionis planetary candidates (6-13 Mj up ) 
in the Palomar and Calar surveys that convey our high cluster 
membership probability criterion and/or fit the near-infrared 
cluster sequence. Another four lack follow-up near-infrared 
photometry. Based on the observed colors, their spectral types 
must range from early-L to mid-L types. From the local 
space number density of L1-L4 dwarfs (Cruz et al. 2003 1, 
we have calculated that about 2-3 field dwarf interlopers of 
similar spectral classifications and brightness in the interval 
z' =21.3-21.9 mag may be contaminating our cr Orionis sur- 
vey. Hence, the relatively high number of L-type cluster 
candidates found in our search can be understood in terms of 
the larger population of cr Orionis as compared to the field. 

3.2. Praesepe and cr Orionis substellar mass functions 

We have obtained the Praesepe and cr Orionis mass functions 
by considering the recent evolutionary models of the Lyon 
group (Baraffe et al. 2003 and references therein), cluster ages 
of 500 and 3 Myr, and distances of 180 and 352 pc, respec- 
tively. For each mass, the theoretical surface temperature and 
luminosity were converted into observables in the manner de- 
scribed in preceeding Sections. Only known and new candi- 
dates conveying our photometric criterion of high-probability 
cluster membership, and candidates that do follow the near- 
infrared photometric sequences of the clusters were taken into 
account to derive the mass functions depicted in Fig. [8] Thus, 
previously known objects like S Ori 42, 44, and 67 are not in- 
cluded in the cr Orionis mass function calculation. Our data 
are plotted as filled dots and histograms, and cover the fol- 
lowing mass intervals: 0.045-0.13 M (Praesepe), and 0.007- 
0.072 M (cr Orionis). We note that all histogram bins of Fig. [8] 
lie within the z'-band completeness of our search and contain 
at least 5 cluster member candidates with the only exception 
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Fig. 8. Praesepe (500 Myr, left panel) and cr Orionis (3 Myr, right panel) mass functions from stars to substellar objects. The 
results from our surveys are depicted with filled dots and histograms. Error bars are Poissonian. In the left panel, we have 
also plotted the Praesepe stellar mass functions from Adams et al. (2002 solid line) and Hambly et al. (Q995a dotted line). In 
both panels, and for comparison purposes, the cr Orionis mass function of Caballero (2006) is shown with a dashed line. It is 
conveniently scaled to the stellar (Praesepe) and substellar (cr Orionis) regimes. Vertical dotted lines correspond to the hydrogen 
(0.075 M ) and deuterium (0.012 M ) burning-mass limits. 



of the most massive bins, which remain unfilled due to non- 
linearity and saturation effects of the CCDs, and the substellar 
bin of Praesepe. 

As discussed by Bejar et al. (2001), the number of 
cr Orionis members appears to increase toward lower masses 
with a smooth transition throughout the deuterium burning- 
mass borderline. Overplotted in Fig. [8] is the cluster stellar 
and substellar mass function recently obtained by Caballero 
( 2006 1. We remark the very good agreement between Caballero 
( 2006 1 and our derivation for the substellar regime. In the com- 
mon mathematical linear representation of the mass function, 
dN '/dm ~ mT a , our data can be fit with an exponent a = 0.6 
where the uncertainty accounts for the different theoretical 
mass-luminosity relationships corresponding to an age interval 
1-10 Myr. 

The substellar census of Praesepe presents a different pic- 
ture. In the following discussion we assume the age of Praesepe 
is ~500Myr. In the left panel of Fig. [3] we have included the 
cluster stellar mass functions from 0.1 up to 1 M obtained 
by Adams et al. (2002i and Hambly et al. H1995at . The for- 
mer authors determined that the great majority of the Praesepe 
stars are mainly contained within a radius of 3.8 deg, while 
the core radius of the cluster is about 1 . 1 deg as derived from 
the King model fitting of the surface density of stars (King 
I1962> . Our survey explores regions within the cluster core ra- 
dius, where we expect the largest number of brown dwarfs to be 
found, and clearly extends the mass functions of Adams et al. 



( 2002) and Hambly et al. ( 1995a) toward the substellar domain. 
For a proper comparison between our Praesepe mass function 
and the results drawn from the literature, we have normalized 
the mass functions of Adams et al. (2002) and Hambly et al. 
( 11995 a) by taking into account the ratio between the total clus- 
ter area and our surveyed area, and the assumption that both 
brown dwarfs and low-mass stars share a similar spatial distri- 
bution in Praesepe. As can be seen from Fig. [8] Adams et al. 
( 2002 ) found that the number of Praesepe stars rises from 1 .0 to 
0.4 M Q , where the mass function distinctly becomes flat down 
to 0.1 M . This contrasts with the continuing rising function 
below 0.4 M a of Hambly et al. dl995at . The two stellar bins 
derived from our search are consistent with the findings of the 
latter authors. 

In the substellar regime, Chappelle et al. ( 2005 ) argued that 
there is a paucity of brown dwarfs in Praesepe leading to a 
turnover of the mass function near the substellar boundary. We 
can test this via the direct comparison of the results of our sur- 
veys in Praesepe and cr Orionis. In our search, there are 5 high- 
probable Praesepe candidates with masses around the substel- 
lar limit (70-90 Mj up ), while only one high-probable candidate 
has appeared at much fainter magnitudes and with a likely mass 
of 52 Mj up at 500 Myr. On the contrary, we estimate that 4-6 
Praesepe brown dwarfs should have been found in the mass in- 
terval 45-65 Mj up if the Praesepe and cr Orionis clusters have 
similar mass functions. Furthermore, the cr Orionis photomet- 
ric sequence of our survey extends well below the deuterium 
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burning mass limit down to about 6 Mj up . The apparent drop 
of a factor of 5 in the number of Praesepe cluster members be- 
tween the substellar limit and the brown dwarf regime is indeed 
indicative of the deficit of massive brown dwarfs in the cen- 
tral regions of Praesepe, confirming the previous conclusions 
of Chappelle et al. ( 2005 ). This result is also in agreement with 
the findings in the Hyades cluster by Dobbie et al. ( 2002 ). 

It also points to a marked discrepancy between the 
<x Orionis and Praesepe substellar mass functions. As can be 
seen from Fig. |8j there is a shallow rising transition between 
stars and brown dwarfs in cr Orionis, while the substellar bin 
(~55 Mj up ) of Praesepe, which is affected by a strong uncer- 
tainty since it is based on one single candidate, appears sig- 
nificantly lower than the contiguous, more massive stellar bin. 
This suggests a sharp change of the Praesepe mass function at 
the substellar boundary to a rapidly decreasing trend for smaller 
masses. The population of brown dwarfs as compared to stars 
seems to be more abundant in cr Orionis than in the central re- 
gions of the Praesepe cluster. 

There are several explanations to account for the discrep- 
ancy between the cr Orionis and Praesepe substellar mass func- 
tions that we will enumerate below. First, our survey covers 
a quite small area of the Praesepe cluster, i.e., it might yield 
results that are non-representative of the entire cluster popula- 
tion. However, we note that other surveys mapping larger areas 
(e.g., Chappelle et al. 2005 1 reach similar conclusions to ours. 
Second, Praesepe is a much older star cluster when compared 
to cr Orionis. Thus, it has likely experienced some mass segre- 
gation and dynamical evolution: the smallest cluster members 
increase their kinetic energy and may have been dispersed from 
the central regions of Praesepe due to encounters and gravita- 
tional interactions with more massive cluster stars. Actually, 
Adams et al. (2002i concluded that there has been some evap- 
oration in Praesepe, i.e., the core of the cluster is depleted of 
low-mass stars. Because of these effects, future searches for 
the Praesepe brown dwarf population should be conducted in 
the halo of the cluster. 

Third, it might be possible that substellar objects cool down 
much faster than predicted by current evolutionary models. 
This would make brown dwarfs fainter than expected at the 
age of Praesepe, i.e., there is a steeper mass-luminosity rela- 
tion between the cluster smallest stars in thermonuclear equi- 
librium and the most massive brown dwarfs. Member candi- 
dates like PraeJ084 11 1.9+192043 would have larger masses. 
The substellar borderline is accurately determined to happen at 
M6-M7 spectral types (-2800 K) in the Pleiades (Stauffer et 
al. H9981 Martin et al. fT998l . which has an age of 120 Myr. 
It is expected that the star-brown dwarf frontier lies at cooler 
types for older ages. Lithium is used as a diagnostic for substel- 
larity (Rebolo et al. H992i : the detection of this atomic feature 
in the spectra of old, ultracool dwarfs implies lithium preser- 
vation and a mass below the substellar limit (see Basri 2000 
for a comprehensive review). As discussed by Kirkpatrick et 
al. l|2000 the presence of lithium in the optical spectra of field 
LO-type dwarfs is ~0%, and slowly rises to 20% by L2, and to 
50-70% by L6. Thus, about 50-70% of the field mid-L type 
dwarfs are bona-fide brown dwarfs for which models predict a 
mass smaller than about 60 Mj up . This suggests that the sub- 



stellar borderline is located between spectral types M6 and L6 
in both the field and the Praesepe cluster. 

Other alternative possibility is the non-universality of the 
initial mass function. However, to address this issue a wide va- 
riety of much younger star clusters have to be studied. On the 
one hand, the Praesepe stellar mass function between 1 and 
0.4 M Q resembles that of the cr Orionis cluster (see Fig. [HJ. On 
the other hand, evidence for mass segregation in the Praesepe 
cluster has been reported in the literature (e.g., Raboud & 
Mermilliod 1998 1, which makes this cluster a non-ideal site to 
investigate the initial mass function. 

Finally, it is possible that the age of the Praesepe cluster is 
significantly larger than 500 Myr. This has an impact on the 
substellar mass function calculation since theoretical brown 
dwarf mass-luminosity relationships depend strongly on age. 
An older age of 800 Myr or higher would bring the substellar 
mass functions of Praesepe and cr Orionis to a better agree- 
ment. We estimate the mass of our coolest Praesepe candidate, 
Prae J0841 1 1.9+192043, to be 57 M Jup at the age of 800 Myr, 
and 63 Mj up at 1 Gyr. However, the mass bin at the substel- 
lar limit roughly changes from 500 Myr to 1 Gyr, suggesting a 
shallow transition from Praesepe very low-mass stars to mas- 
sive brown dwarfs, as in cr Orionis. The various determinations 
of the age of Praesepe available in the literature are generally 
in the range 500-900 Myr, and recent determinations based on 
the isochrone-fitting procedure yield a likely age of 800 Myr 
(Kharchenko et al. 2005 1. For this and older ages, the bulk of 
the cluster brown dwarf population lies at magnitudes fainter 
than those of Prae J0841 1 1.9+192043 and will show colors 
typical of late-L and T-types. Future searches should take this 
possibility into account. 

4. Summary 

We have conducted deep photometric optical surveys in the 
Praesepe (0.5-1 Gyr) and cr Orionis (~3 Myr) clusters us- 
ing the Sloan i' and z' filters, large format multi-CCD cam- 
eras, the 3.5 m telescope on the Calar Alto Observatory and 
the 5-m Hale Telescope on the Palomar Observatory. Our 
main objective is the identification of the reddest population 
(if — z' > 1 mag) of each cluster and place constrains on the bot- 
tom end of the Praesepe and cr Orionis mass functions. A total 
area of 1177 (Praesepe) and 1122 (cr Orionis) arcmin 2 were 
explored in the cluster central regions down to limiting mag- 
nitudes of i' =24.5 and z' =23.5 mag (completeness is about 
one magnitude brighter). Our search is sensitive to Praesepe 
and cr Orionis members with masses in the intervals 0.05- 
0.13 M Q and 0.006-0.072 M Q , respectively. Follow-up near- 
infrared JHK S photometry has allowed us to discuss the cluster 
membership of various candidates identified from the optical 
surveys. 

In cr Orionis, three new planetary-mass candidates 
(SOriJ053844.5-025514, S Ori J053932.4-025220, and 
SOriJ054008. 5 -024551) with estimated masses (0.006- 
0.010 M Q ) below the deuterium burning-mass limit have been 
discovered. They fit the optical and near-infrared sequence 
of the cluster, suggesting that they are likely members of 
cr Orionis. From our survey, we derive a rising substellar mass 
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function which extends from the substellar borderline to the 
planetary regime (a = 0.6 AN /AM ~ M~ a ), in agreement 
with our previous derivations (Bejar et al. 120011 Caballero 
120061 . 

In the Praesepe cluster, we have found one red candidate 
(Prae J0841 1 1.9+192043) whose optical and near-infrared col- 
ors coincide with our expectations for true substellar members. 
Its mass and spectral type are estimated at 50-60 M Sup and L4- 
L6, respectively. It is the coolest and least massive Praesepe 
candidate reported to date. At the youngest possible ages of 
Praesepe (500-700 Myr), our analysis of the cluster mass func- 
tion suggests a sharp change at the substellar boundary to a 
rapidly decreasing trend for brown dwarf masses, which clearly 
contrasts with the <x Orionis mass function. This trend dimin- 
ishes to a shallow transition between very low-mass stars and 
massive brown dwarfs when using steeper mass-luminosity re- 
lationships like those given by the models of 800 Myr and 
higher. 
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Fig. 9. Finder charts (z'-band, 2' x 2' in extent) of Praesepe candidates labeled with a "high" photometric cluster membership 
probability. [Provided in jpeg format in astro-ph.] 

Fig. 10. Finder charts (z'-band, 2' x 2' in extent) of Praesepe candidates labeled with a "low" photometric cluster membership 
probability. [Provided in jpeg format in astro-ph.] 

Fig. 11. Finder charts (z'-band, 2' x 2' in extent) of cr Orionis candidates labeled with a "high" photometric cluster membership 
probability. [Provided in jpeg format in astro-ph.] 

Fig. 12. Finder charts (z'-band, 2' x 2' in extent) of cr Orionis candidates labeled with a "low" photometric cluster membership 
probability. [Provided in jpeg format in astro-ph.] 
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